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ABSTRACT 

We present multi-epoch non-redundant masking observations of the T Cha transition disk, taken at 
the VLT and Magellan in H, Ks, and L' bands. T Cha is one of a small number of transition disks 
that host companion candidates discovered by high-resolution imaging techniques, with a putative 
companion at a position angle of 78°, separation of 62 mas, and contrast of A L' = 5.1 mag. We 
find comparable binary parameters in our re-reduction of the initial detection images, and similar 
parameters in the 2011 L', 2013 NaCo L', and 2013 NaCo Ks data sets. We find a close-in companion 
signal in the 2012 NaCo L' dataset that cannot be explained by orbital motion, and a non-detection in 
the 2013 MagAO/Clio2 L' data. However, Monte-carlo simulations show that the best fits to the 2012 
NaCo and 2013 MagAO/Clio2 followup data may be consistent with noise. There is also a significant 
probability of false non-detections in both of these data sets. We discuss physical scenarios that could 
cause the best fits, and argue that previous companion and scattering explanations are inconsistent 
with the results of the much larger dataset presented here. 


1. INTRODUCTION 

Since their discovery, transition disks have been re¬ 
garded as natural laboratories for the study of proto¬ 
planetary disk evolution and perhaps planet formation. 
These objects’ spectral energy distributions (SEDs) lack 
near- t o mid-infrared emissio n, yet display a far-infrared 
excess (|Strom et al. ||1989|). Initial studies attributed 
these SED features to a lack of warm dust at inner, 
AU-scale radii, suggesting that they were “in transition” 
from protoplanetary disks with excess throughout the in¬ 
frared, to de bris disks with only very weak far-infrared 
excess (e.g., Lin fe Papaloizou 1986 1993 |Bryden et 


al. 


1999 


Calvet et al. _ 5nU2|). More recenf modeling 


of ilpitzer spectra (e.g., [Calvet et al. || 20051 |Brown et 
al. 2007 Espaillat et al. 2007a|b"Merin et al. 2010|) 
also associated the mid-infrared deficits with disk cavities 
or ga ps on AU scales. Followup submillimeter imaging 
’e.g., B rown et al.~||2009 Hughes et al. |[2009 Andrews 


et al. |2Ull[) has directly confirmed the presence of these 


features. 

Studies have identified several processes that could 
play a role in the formation of gaps and cavities, includ¬ 
ing photo-evaporative winds, grain growth, and dynam¬ 
ical interactions with companions. While photoevapora- 


in the 

inner disk, ( 

Clarke et al. 

2001 

Alexander et al. 

2006 

), accounting tor radial dri 

t ol solids can lead to 


dissipation of du st at small radii as well (Alexander & 


Armitage 2007). Furthermore, X-ray winds may drive 
disk depletion at a faster rate than UV winds, suggesting 
that X-ray photoevaporation could clear our inner disk 
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radii more efficiently (e.g., 

Ercolano et al. 

2008 

Drake 

et al. 

2009 

Owen et al. 

2010). However, t 

ie cav- 


are too large compared to their X-ray luminosi ties to be 


consistent wi t h clearing by photoe vaporation (Andrews 


et al. 2011 Owen et al. 


2011). Furthermore, mea¬ 


surements of outer disk masses are too large compared 
to results of simulation s that cause disks to go through 


an “inner hole” phase (Alexander et al. 2006). While 


photoevaporation could explain some inner clearings, it 
alone cannot be responsible for transition disk structure. 

Rather than clearing away disk material to lower the 
infrared emission, grain growth decreases its emissivity 


(Draine et al. 2006 D’ Alessio et al. 


grains to mm/cm sizes can create SEP defici ts compara- 


2006). Growing 


ble to those observed in transition disks (e.g.JDullemon d 
fe Dominik ||2005 Tanaka et al. 2005| |Birnstiel et al. 


2012|). However, disk evolution simulations by Birnstiei 


et al. (2012) failed to generate the particle size distribu¬ 
tions required to produce mm wavelength cavities. Ad¬ 
ditionally, soon after growing from ~ 1/rm size to mm 
size, direct collisions b etween silicate particles could be¬ 
come destructive (e.g., Windmark et al. ] |2012 1, and the 
resulting smaller particles could then produce emission 
to fill in the SED deficit. This suggests that, while grain 
growth must impact disk evolution at some level, this 
process alone cannot shape transition disk cavities. 

Dynamical interactions with companions are the best 
explanation to date for forming disk gaps. Models have 
demonstrated that stella r mass companions can ope n 


cavities in disks (e.g., Artymowicz & Lubow 1994 


Pichardo et al. 20081, and s ome observed disk gaps 
such as th ose in CoKu Tau 4 dlreland fe Kraus ||2008|) 

Uchid; 


HD988 00 (|Furlan et al. 2007), and Hen 3-600 


et al. |[2004), have been associated with stellar mass bi¬ 


naries. However, high resolution imaging has ruled out 
companions with masses higher than ~ 20 - 30 M j,,,, for 


approximately half of the know n transition disks (Kraus 


et al. |20TT Evans et al. 2012J). This leaves the exciting 
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possibility that planetary mass companions are clearing 
out cavities and gaps, accreting mat erial that would hav e 
otherwise fallen onto the star (e.g., Najita et al. J2007). 


Simulations have shown that tidal interactions with a 


(e.g., 

Lin & Papaloizou |1986 

Bryden et al. 

1999 

Crida 

et al. 12006). 


shaped by a substellar mass companion, T Chamaeleon- 
tis (T Cha). T Clra is a G8 type, 1.5 M 0 star, first 
categorized as a weak-line T Tau ri star due to its low 


Ha equivalent width of < 10 A (Alcala et al. 1993). 
This suggested it had entered the final stages of accre¬ 
tion. The classification conflicted with T Cha’s infrared 
excess, thought to result from its outer disk; this spec¬ 
tral feature would place it as a classical T Tauri star. 
Later observations showed that it is in fact a transition 
disk object, perhaps in the intermediate stages between 
a protop lanetary dis k and a disk-free planetary system. 


Brown et al. (2007) found that T Cha’s SED could be 
reproduced by a disk with a gap between 0.2 and 15 AU, 


and SED modeling by Olofsson et al. (2011) supported 
this, with a best fit gap extending from 0.17 to 7.5 AU. 

Imaging observations suggested the presence of a com¬ 
panion of L' contrast 5.1 mag at a se paratio n of 62 m as 
- 6.7 AU at a distance of 108 pc - jHu elamo et al. 
2011|), within T Cha’s disk gap. However, Ulolsson et af. 


11 

al: 


2013) showed that the detected signal could be modeled 


almost equally well by an asymmetry caused by forward 
scattering from the upper layers of the outer disk. Ob¬ 
serving orbital motion of the companion candidate would 
distinguish between these two scenarios. To this end, we 
acquired new observations of T Cha with the Magellan 
AO (MagAO) system. We also present a re-analysis of 
the original discovery data from VLT/NaCo, as well as a 
new analysis of previously unpublished NaCo data from 
the Very Large Telescope (VLT) archive. 


2. EXPERIMENTAL SETUP 

Non-redundant masking (NRM) transforms a filled 
aperture into a sparse interferometric array using a pupil- 
plane mask. While blocking the majority of the light, this 
provides much better knowledge of the PSF than a con¬ 
ventional telescope. A resulting image (called an interfer- 
ogram) then shows the interference fringes formed by the 
mask, and subsequent image reconstruction or model fit¬ 
ting relies on quantities calculated from its Fourier trans¬ 
form. Since the mask is non-redundant (no two baselines 
have the same length and orientation), each baseline has 
unique (it, v) coordinates. The symmetry of the Fourier 
transform means that identical information for a baseline 
can be found in two points in (it, v ) space - at (it, v) and 
(— u, —v). The finite size of the mask holes, as well as the 
width of the filter bandpass, causes this information to 
spread out. This means that the Fourier transform will 
have several distinct “splodges,” two coming from each 
mask baseline. 

Using the Fourier transform at the locations of these 
splodges, we find the complex visibility for each base¬ 
line, which has the form Ae where A is the ampli¬ 
tude and (/> the phase. Since atmospheric and instru¬ 
mental effects corrupt the complex visibilities, we calcu¬ 
late two other quantities, squared visibilities and closure 
phases. Squared visibilities measure the power in the 


Fourier transform as a function of baseline length, and 
closure phases are sums of phases around three baselines 
that form a triangle. Closure phases eliminate atmo¬ 
spheric and instrumental phase offsets that corrupt mea¬ 
surements taken using single baselines. These obey the 
relation: 


<F cp = c/)(ui,v 1 ) + (j)(u 2 ,v 2 ) +4 '(u 3 ,v 3 ) (1) 

where Ui and Vi are the sampling coordinates of the 
i th baseline in the Fourier plane. An N-hole mask will 
provide (^) baselines and visibilities, and (^) closure 
phases, ( A ^" 1 ) of which are independent. 


3. OBSERVATIONS 

3.1. New Magellan/MagAO/Clio2 Data 

We observed T Cha and two unresolved calibrators 
HD101251 and HD102260 using the 6.5 m Clay tele¬ 
scope , MagAO adap tive optics system (Close et al. 2013 

- 


MorzinskiAtaL|26l4 ), and Clio2 science camera 


breed 


et af. 2004 " bivananc am et al. "|2006| on 2013 ApHTST 
A six-hole non-redundant mask was mounted in a pupil 
plane filter wheel in Clio2. Table [T] lists the parameters 
of the Clio2 mask. We used two calibrators to lessen 
the probability that detected signals were being injected 
by a resolved or binary calibrator. Our exposure times 
for T Cha, HD101251, and HD102260 were 1.0 s, 0.9 
s, and 1.3 s, respectively. We took 4 to 6 50-frame data 
cubes (called “visits” in Table[2]) before switching objects 
according to the pattern target-call-target-cal2. We ac¬ 
quired 6 visits for T Cha, 2 visits for HD101251, and 3 
visits for HD102260. These resulted in 28.3 minutes, 9 
minutes, and 15.16 minutes of total integration, respec¬ 
tively. We observed in L', A c = 3.78 pan, and the total 
change in sky rotation angle was 47°, shown in Figure [T] 


3.2. Previously Published 2010 VLT/NaCo Data 
3.2.1. 2010 March 17 Ll Observations 


Th e in itial detection of the T Cha companion (Hue- 
lamo et al . 12011 ) resulted from NRM observations taken 


at th e V LT using NaCo (Penzen 
2003). The parameters 


1993 Rousset et al. 
or the iSlaCo mask are fisted 
in Table |T| In order to verify our reduction pipeline, 
we downloaded and re-reduced the previously published 
observations from the archive. This dataset included ob¬ 
servations of T Cha and a single, unresolved calibrator, 
HD102260, with 9 visits to T Cha and 10 to HD102260. 
The observations followed the pattern ...cal-target-cal..., 
dithering so that each 100-frame data cube placed the 
interferogram on one of the detector’s quadrants. These 
data resulted in 48 minutes of integration on T Cha and 
53.33 minutes on HD102260. Images were taken in L', 
A c = 3.8 /im, and the total change in sky rotation angle 
was 62°, shown in Figure [l] 


3.2.2. 2010 July 1 Ks Observations 

The detection of T Cha in L' was accompanied by a 
non-detection in Ks (A c = 2.18/im), from data taken 
at the VLT using NaCo in July 2010. We reduced this 
archival dataset, which includes observations of T Cha, 
HD102260, and HD101251. Three visits were made to 
T Cha, totaling 20 minutes of integration. For each cal¬ 
ibrator, 1 visit consisting of 6.7 minutes of integration 
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TABLE 1 
Mask Parameters 


Mask 

Wholes 

N b 

N cp 

Nk p 

Baseline Range (m) 

Throughput (%) 

NaCo 

7 

21 

35 

15 

1.77 - 6.45 

16% 

MagAO / Clio2 

6 

15 

20 

10 

1.68 - 5.02 

11% 


NaCo L' 2010 



NaCo L' 2011 NaCo L' 2012 



MagAO/Clio2 L' 2013 NaCo Ks 2013 



u (m) 


Fig. 1.— Sky rotation comparison in (u,v) space for all L' data sets and 2013 NaCo Ks data. The NaCo 2010 L' data had a change in 
sky rotation of ~ 62°. The total sky rotation in the NaCo 2011 L' data was ~ 39°. The NaCo 2012 L' data had ~ 19° of sky rotation, and 
the NaCo 2013 L' data had ~ 56°. The MagAO/Clio2 2013 data had a change of ~ 47°. Due to its smaller aperture, the MagAO/Clio2 
baselines are shorter than for NaCo. Lastly, the NaCo 2013 Ks data had ~ 71° of sky rotation. 


was made. The change in sky rotation angle for these 
data was 19°, and the dithering pattern was identical to 
the 2010 L ' NaCo data. These data were too noisy to de¬ 
tect companion signals comparable to the Huelamo et al. 
I (2011) binary, so we include their discussion in Appendix 

3.3. Unpublished VLT/NaCo Data 

While searching for published data to verify our 
pipeline, we found additional observations taken using 
VLT/NaCo from 2011 - 2013. These include L' and Ks 
observations from March 14-15, 2011, L' observations 
from March 8, 2012, and L/, Ks, and H band observa¬ 
tions from March 25-27, 2013. We include a description 
of each dataset in this section. However, the scatter in 
the 2011 Ks and 2013 H band observations would wash 
out companion signals of interest. For this reason, we 


include only a short discussion of the results from these 
data in Appendix [A] 

3.3.1. 2011 March 14 L' Observations 

We reduced archival L ' VLT/NaCo data taken in 2011, 
which included observations of T Cha and the same cal¬ 
ibrator as the 2010 L' dataset, HD102260. These data 
consisted of seven visits to T Cha, and nine visits to 
HD102260, resulting in 37.33 minutes of integration on 
T Cha and 48 minutes on HD102260. The dithering pat¬ 
tern was the same as for the 2010 L' NaCo data, and 
the total change in sky rotation angle was 39°, shown in 
Figure []] 

3.3.2. 2011 March 15 Ks Observations 

These archival data include observations of T Cha and 
HD102260 taken in Ks. The dithering pattern was iden- 
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TABLE 2 

Summary of Observations 


Target 

Right Ascension 
(hh mm ss.sss) 

Declination 
(dd mm ss.sss) 

t int N frames a b Total Time Seeing r 0 

(s) (min) (arcsec) (ms) 



L' 

Observations 


2010 Mar 17: VLT/NaCo 

T Cha 11 57 13.550 

HD102260 11 45 13.822 

-79 21 31.537 
-78 36 58.633 

0.4 

0.4 

800 

800 

9 

10 

48 

53.3 

0.6 

8 

2011 Mar 

14: VLT/NaCo 








T Cha 

11 57 13.550 

-79 21 31.537 

0.4 

800 

7 

37.3 

2.0 

1.5 

HD102260 

11 45 13.822 

-78 36 58.633 

0.4 

800 

9 

48 



2012 Mar 8: VLT/NaCo 








T Cha 

11 57 13.550 

-79 21 31.537 

0.4 

800 

5 

26.7 

1.0 

6 

HD101251 

11 37 49.220 

-79 14 31.604 

0.4 

800 

2 

10.7 



HD102260 

11 45 13.822 

-78 36 58.633 

0.4 

800 

2 

10.7 



2013 Mar 25: VLT/NaCo 








T Cha 

11 57 13.550 

-79 21 31.537 

0.4 

800 

2 

10.7 

0.75 

7 

T Cha 

11 57 13.550 

-79 21 31.537 

0.3 

1057 

7 

37.0 



T Cha 

11 57 13.550 

-79 21 31.537 

0.2 

1405 

1 

4.7 



HD102260 

11 45 13.822 

-78 36 58.633 

0.4 

800 

1 

5.3 



HD102260 

11 45 13.822 

-78 36 58.633 

0.3 

1057 

4 

21.1 



HD102260 

11 45 13.822 

-78 36 58.633 

0.5 

1057 

1 

8.8 



HD101251 

11 37 49.220 

-79 14 31.604 

0.4 

800 

1 

5.3 



HD101251 

11 37 49.220 

-79 14 31.604 

0.3 

1058 

2 

10.6 



HD101251 

11 37 49.220 

-79 14 31.604 

0.15 

1687 

1 

4.2 



2013 Apr 5: Magellan/MagAO/Clio2 







T Cha 

11 57 13.550 

-79 21 31.537 

1.0 

300 

5 

28.3 

0.6 

N/A 

T Cha 

11 57 13.550 

-79 21 31.537 

1.0 

200 

1 

28.3 



HD101251 

11 37 49.220 

-79 14 31.604 

0.9 

300 

2 

9 



HD102260 

11 45 13.822 

-78 36 58.633 

1.3 

300 

1 

15.16 



HD102260 

11 45 13.822 

-78 36 58.633 

1.3 

200 

2 

15.16 




Ks Observations 


2010 Jul 1: 

T Cha 

HD101251 

HD102260 

VLT/NaCo 

11 57 13.550 

11 37 49.220 

11 45 13.822 

-79 21 31.537 
-79 14 31.604 
-78 36 58.633 

0.5 

0.5 

0.5 

800 

800 

800 

3 

1 

1 

20 

6.7 

6.7 

1 

4 

2011 Mar 15: VLT/NaCo 








T Cha 

11 57 13.550 

-79 21 31.537 

2.0 

200 

6 

40 

0.75 

3 

HD102260 

11 45 13.822 

-78 36 58.633 

2.0 

200 

6 

40 



2013 Mar 26: VLT/NaCo 








T Cha 

11 57 13.550 

-79 21 31.537 

1.0 

280 

2 

9.3 

1 

4 

T Cha 

11 57 13.550 

-79 21 31.537 

0.8 

427 

4 

22.8 



T Cha 

11 57 13.550 

-79 21 31.537 

0.7 

497 

3 

17.4 



T Cha 

11 57 13.550 

-79 21 31.537 

0.5 

807 

2 

13.5 



HD102260 

11 45 13.822 

-78 36 58.633 

1.0 

350 

1 

5.8 



HD102260 

11 45 13.822 

-78 36 58.633 

0.5 

707 

2 

11.8 



HD102260 

11 45 13.822 

-78 36 58.633 

0.4 

707 

1 

4.7 



HD102260 

11 45 13.822 

-78 36 58.633 

0.3 

1057 

1 

5.3 



HD101251 

11 37 49.220 

-79 14 31.604 

0.2 

1407 

3 

14.1 



HD101251 

11 37 49.220 

-79 14 31.604 

0.11 

2782 

1 

5.1 




H Observations 


2013 Mar 27: VLT/NaCo 


T Cha 

11 57 13.550 

-79 21 31.537 

0.4 

605 

1 

4.03 

T Cha 

11 57 13.550 

-79 21 31.537 

0.4 

847 

1 

5.64 

T Cha 

11 57 13.550 

-79 21 31.537 

0.5 

707 

4 

23.6 

HD102260 

11 45 13.822 

-78 36 58.633 

0.15 

2107 

3 

15.8 

HD101251 

11 37 49.220 

-79 14 31.604 

0.11 

3157 

2 

11.6 

HD101251 

11 37 49.220 

-79 14 31.604 

0.11 

1353 

1 

2.5 


a Number of frames in each visit 

b Each visit consists of all images taken before switching between target and calibrator. 
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tical to the 2010 L' NaCo data. Six visits were made to 
each object, resulting in 40 minutes on target and cali¬ 
brator. The total change in sky rotation for these data 
was 38°, shown in Appendix [A] 

3.3.3. 2012 March 8 L' Observations 

We reduced archival VLT/NaCo L' data taken on 
March 8, 2012. These include observations of T Cha, 
HD102260, and HD101251. Five visits were made to T 
Cha and 2 to each calibrator. This resulted in a total of 
26.7 minutes of integration for T Cha, and 10.7 minutes 
for each calibrator. The visits were dithered so that the 
image fell on a different detector quadrant during neigh¬ 
boring sets of 100 exposures. The total change in sky 
rotation for this dataset was 19°, shown in Figure [T] 

3.3.4. 2013 March 25 L' Observations 

The 2013 archival L' VLT/NaCo data include obser¬ 
vations of T Cha, HD101251, and HD102260. Ten visits 
were made to T Cha, with 52.3 total minutes of inte¬ 
gration. A total of 4 visits were made to HD101251, 
resulting in 20.1 minutes of integration, and the 6 visits 
to HD102260 yielded 35.3 minutes of integration time. 
These data had a change in sky rotation of 56°, shown in 
Figure [l] Table [2] details the individual visits to all three 
objects, which did not have identical exposure times. 

3.3.5. 2013 March 26 Ks Observations 

We also present archival VLT/NaCo data taken in Ks 
band. The 11 visits to T Cha resulted in a total of 63 
minutes of integration. HD101251’s 4 visits yielded 19.2 
minutes of total exposure time, and the total amount of 
integration for the 5 visits made to HD102260 was 26 
minutes. The change in sky rotation for this dataset was 
71°, shown in Figure [T| See Table [2] for details of the 
individual visits to each object. 


3.3.6. 2013 March 27 H Observations 

The last archival dataset was taken on March 27, 2013 
in H band (A c = 1.65/im). The total integration time 
for T Cha’s 6 visits was 33.3 minutes. For HD101251, 2 
visits were made, resulting in 14.1 minutes of integration. 
Lastly, the 3 visits to HD102260 total 15.8 minutes of 
integration time. Table [2] details the individual visits 
made to each object. The change in sky rotation for 
these observations was 39°, shown in Appendix [A] 


4. DATA REDUCTION 

We have developed a suite of software in Python to per¬ 
form basic data reduction, calibration, and visibility and 
closure phase calculations. We first flat-field and bad- 
pixel correct all images. For a given set of two dithers, 
we perform sky subtraction for one position by taking the 
median of all images in the other dither position and sub¬ 
tracting the median from each image. We then apply a 
super-Gaussian window function to reduce the noise asso¬ 
ciated with the low-signal edges of the interferogram and 
spatially filter the data. A super-Gaussian has the form 
exp(— kx 4 ); we choose k such that the half width at half 
max is \/d S ub± wh ere d su h is the mask sub-aperture di¬ 
ameter (e.g., Bernat |2 012). After windowing, we Fourier 
transform the data. 



-10 -5 0 5 10 

Kernel Phase (°) 


Fig. 2.— Normalized histograms of uncalibrated kernel phases 
for L/ and 2013 Ks data. For a subset of each dither (chosen so that 
equal amounts of integration came from all observations), we sub¬ 
tract the mean kernel phase from each individual measurement to 
generate the histograms shown above. The snapshot kernel phase 
errors are much lower in the MagAO/Clio2 data than in the NaCo 
data sets, indicating lower levels of random noise. 


Next, we calculate squared visibilities. We first sim¬ 
ulate an interferogram and resulting Fourier transform 
using the locations and sizes of the holes in the mask, 
along with the observation’s wavelength. This provides 
us with the pixel locations of the splodges (see § [2]) in the 
Fourier transform. We then square the Fourier transform 
of the data, and sum all pixels within the splodges corre¬ 
sponding to each baseline, normalizing by the total power 
in the interferogram. 

The typical uncertainty due to random errors for the 
visibilities ranges between 0.02 and 0.12 for all L' and 
2013 Ks data sets. For comparison, a binary with 
separation 62 mas and contrast A L' = 5.1 mag pro¬ 
duces a change in visibility between the shortest and 
longest baseline of approximately 0.035 for both the Ma- 
gAO/Clio2 and NaCo masks. While some of the followup 
data sets’ random visibility errors are smaller than this 
signal, visibilities’ dependence on the AO system (e.g., 
. |2011 Kraus & Ireland |2012) renders them 
harder to calibrate. Differences in AO performance over 
the night, or between target and calibrator observations 
can introduce additional error. To investigate this, we di¬ 
vide the visibilities for each set of two adjacent calibrator 
scans; the calibrated visibilities for a point source should 
be equal to 1 for all baselines. We then take the scatter in 
these calibrated visibilities as an estimate for the system¬ 
atic uncertainties in the target visibilities. This results in 
calibrated visibilities with scatters ranging between 0.04 
and 0.11. Due to th ese large uncertainties, as in previous 


NRM stud ies (e.g., Huelamo et al. 2011 Kraus & Ire- 
land |[20l2) we restrict our binary fitting to phases only. 


rather than including the squared visibilities. 

For each triangle, we form the bispectrum by multiply- 
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Kernel Phase (°) 


Fig. 3. — Normalized histograms of each set of calibrated L' kernel phases as well as NaCo 2013 Ks kernel phases, with their best fit 
Gaussian distributions over plotted. The Gaussian distributions were used to generate noise realizations for the simulations described in 

Section [3 


ing the complex values in the Fourier transform at three 
(zz, v ) coordinates. We then average over the individual 
frames, and take the phase of the average bispectrum to 
be our closure phase. Of the (^) closure phases, only 
( JV 9 ~ 1 ) are independent. For this reason, we perform fits 
on kernel phases, linearly independent combinations of 
closure phases (see Martinache 2010, Kraus & Ireland 
2012, and Ireland 2013) . W e find our kernel phases in a 


way similar to Martinache (2010). Appendix [C| gives a 


detailed description of our projection method. 

Figure [2] shows histograms of the kernel phases for 
individualimages in the Ks and L' data sets, and can 
be taken to represent a comparison of the snapshot ker¬ 
nel phase errors for the different data sets. For a given 
dither, we calculate the ( Ar “ 1 ) mean kernel phases. We 
then calculate the kernel phases for every individual im¬ 
age in the dither and subtract the mean. We use the same 
total integration time, 39 s, to calculate the mean ker¬ 
nel phase. This process yields the distributions shown in 
Figure [2] The snapshot errors in the 2013 MagAO data 
(er = 0.73°) are much smaller than the VLT/NaCo data 
sets (a = 2.85 — 4.38° for L ' kernel phases.) We also 
compare the MagAO/Clio2 and NaCo data with differ¬ 
ent amounts of time-averaging. As long as these aver¬ 
ages are performed within a single dithering sequence, 
the chosen interval does not change the noise levels sig¬ 
nificantly. Due to their high scatter and the resulting 
unreliable fits we show kernel phase histograms for the 
2010 and 2011 Ks as well as the 2013 H band data in 
Appendix [A] 

Random sources of noise associated with both AO per¬ 
formance and observing conditions cause snapshot kernel 


phase errors. For exposures much longer than the inverse 
of the AO system band width, closu re phase errors scale 
in the following way (Ireland ||2013): 


Gcp 



( 2 ) 


where a$ is the phase noise on each sub-aperture in the 
closing triangle, f c the cutoff frequency below which pis¬ 
ton noise is white, and T the exposure time. The kernel 
phase errors will then be the projection of the closure 
phase errors. If we assume er^ and f c are equal for all 
data sets, the ratio of two observations’ closure phase 
errors scales with the ratio of their exposure times. Us¬ 
ing 0.2s - 0.5s for the NaCo data, and 1.0s for the Ma- 
gAO/Clio2 data, the ratio of the MagAO/Clio2 to NaCo 
closure phase errors should range from 0.44 to 0.71. The 
ratio of the MagAO/Clio2 and NaCo L' observed errors 
ranges between 0.27 and 0.18 depending on the NaCo 
L' dataset. Hence, exposure time alone cannot explain 
the discrepancy between the MagAO/Clio2 and NaCo L' 
snapshot errors. 

Other random sources of noise include photon, back¬ 
ground, and read noise, which lead to closure phase errors 
(following |Ireland ||2013 ): 

a ° v ~ Wv^ l '^ Np + Nb + n P a ro) . ( 3 ) 

where 1W is the number of holes in the mask, V the fringe 
visibility, N p the total number of photons, A), the number 
of background photons, n p the number of pixels, and cr ro 
the read noise. If we assume we are in the photon-noise 
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regime, the closure phase error simplifies to 


Nh /T5 

° cv ~ V V 


( 4 ) 


The number of photons is N p = F target x A te i x / mas fc x T, 
where F tar get is the photon flux from the source, A te i the 
telescope collecting area, fmask the fraction of light al¬ 
lowed through by the mask, and T the exposure time. 
The relevant values for the MagAO/Clio2 observations 
are A tel ~ (6.5m) 2 , f mask = 0.11, T = Is, and N h = 6 . 
For the NaCo observations, A te i ~ ( 8 . 2 m) 2 , f mas k = 
0.16, T = 0.2s — 0.4s, and N k = 7. Assuming that the 
fringe visibilities for the two data sets are approximately 
equal, taking the ratio of the MagAO/Clio2 to NaCo clo¬ 
sure phase errors gives 0.58 — 0.92. Thus, the differences 
in telescope and observing parameters cannot fully ex¬ 
plain the lower MagAO/Clio2 snapshot errors. Better 
AO performance by MagAO, which leads to lower val¬ 
ues for 00 and higher values of f c , could be one cause 
for this discrepancy. Additionally, the smaller holes (as 
evinced by the lower throughput in Table IT]) of the Clio 2 
mask mean that the 2013 observations are less redundant 
than the NaCo observations. This could also reduce the 
snapshot errors for the MagAO/Clio2 data. 


5. CALIBRATION 

Due to atmospheric and instrumental systematics, the 
mean kernel phases themselves (subtracted off in Figure 
U|) can vary substantially throughout the observations. 
To take these effects into account, we subtract our un¬ 
resolved calibrator kernel phases from our target kernel 
phases. We do this in several ways. 

To apply a simple nearest-neighbor calibration, we first 
calculate the time between a given target scan and all cal¬ 
ibrator scans (At). We then average all calibrator scans, 
weighting by At” 10 , and subtract the weighted-average 
calibrator from the target scan. To use information from 
all of the calibrator scans, rather than limiting ourselves 
to only the nearest-neighbor, we calculate an average cal¬ 
ibrator, weighting the scans by At” 1 . Finall y, we a pply 
a more opt imized weighting, similar to LOCI (Lafreniere 
et al. |2007 ) techniques in direct imaging data reduction 


and the calibration strategy adopted in Kraus & Ireland 
(2012) and Ireland (2013). 

For the LOCI-like calibration scheme, we find the lin¬ 
ear combination of calibrator scans that minimizes the 
sum of the target’s squared kernel phases. This is equiv¬ 
alent to minimizing the y 2 for the null model. While this 
calibration scheme provides the lowest scatter and thus 
highest signal to noise, it can also subtract signal from 
the measurements. For this reason, it is often applied 
iteratively, minimizing the y 2 for the null model initially 
and then minimizing the y 2 for the best-fit model until 
the best-fit converges. We also LOCI calibrate without 
iteration, minimizing the y 2 of the At -1 model. Both 
LOCI schemes remove signal, and do not always give 
consistent results. Therefore, we focus on the fits to 
simpler, neighbor-like calibrations in the subsequent sec¬ 
tions. The behavior of the LOCI calibration method will 
be detailed in a future paper. 

The histograms for the L' and 2013 Ks data sets are 
shown in Figure [3j While the MagAO/Clio2 data have 
lower snapshot kernel phase errors than the NaCo L ' ob¬ 


servations (see Figure [2]) , the systematics in the three 
data sets are such that the scatter in the mean, cali¬ 
brated kernel phases are quite similar. We speculate that 
the greater temporal spacing between target and calibra¬ 
tor observations in the MagAO/Clio2 observations could 
cause this. The average time between target and cali¬ 
brator scans in the these data is 8.8 minutes. For the 
NaCo observations, the mean time between target and 
calibrator scans for the L' data ranges between 4.2 and 
6.9 minutes. 


5.1. Consistency Checks 

To check whether one calibrator could be contaminat¬ 
ing the MagAO/Clio2 data and increasing the kernel 
phase scatter, we calibrated our T Cha kernel phases 
using each calibrator individually. Our dithering pattern 
alternated calibrator observations between target obser¬ 
vations, resulting in some target scans being closer in 
time to one of the two calibrators. Using only one cali¬ 
brator star increased the scatter in the calibrated kernel 
phases. However, both single-star calibrations yielded 
kernel phases with nearly identical standard deviations 
(approximately 1°). This similarity suggests that neither 
calibrator is contaminating the dataset, and the increase 
in noise compared to the two-star calibration scheme 
highlights the need for calibrator scans taken close in 
time to the target observations. Furthermore, best fits 
to these calibrated data are consistent with those for the 
fully calibrated data. 

For the fully-calibrated MagAO/Clio2 L' data, we com¬ 
pare the scatter in scans taken during the first half of 
the night to the scatter for those taken during the sec¬ 
ond. The kernel phase standard deviations in this test 
are nearly identical (0.66° and 0.68°). This suggests that 
the calibration quality did not change significantly dur¬ 
ing the observations. 

We carried out calibrator tests for the NaCo 2012 L' 
and 2013 L' and Ks data sets, which, unlike the 2010 and 
2011 data, included observations of two calibrator stars. 
We calibrated each calibrator star using the other with a 
At ” 1 method, and then fit a binary model to the kernel 
phases. In all three data sets, neither calibrator star’s 
kernel phases show clear signs of a companion; simula¬ 
tions show that their best fits are consistent with noise. 

We also calibrated the NaCo 2012 and 2013 T Cha 
kernel phases using each calibrator separately, and then 
fit the resulting kernel phases. Using individual calibra¬ 
tors did not change the best fits significantly for any of 
these three data sets. Additionally, the scatters for these 
calibrations were similar. Using only HD101251 and 
HD102260, respectively, the standard deviations were 
0.39° and 0.36° for 2012 L', 0.28° and 0.33° for 2013 
L', and 0.45° and 0.36° for 2013 Ks. 

We also check that the scatters in these three NaCo 
data sets are comparable for both halves of the night. 
The 2012 L' data from the first and second halves of 
the observations have scatters of 0.37° and 0.36°, respec¬ 
tively. For the 2013 L' observations, the kernel phase 
standard deviation for the first half of observations is 
0.32°, while for the second half it is 0.27°. The change 
throughout the night for the 2013 Ks observations was 
slightly larger, with a standard deviation of 0.45° for the 
first half, and 0.31° for the second half. 

Lastly, we compare the scatter in the calibrated kernel 
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TABLE 3 

Binary Grid Parameter Space 


Parameter 

Minimum 

Maximum 

Step Size 

P.A. (°) 

-180 

180 

1 

s (mas) 

0 

700 

5 

A (mag) 

3 

7 

0.05 


phases (Figure pi) with changes in average seeing and co¬ 
herence time (Table [2]) between observations. In general, 
observations with longer coherence times and lower wind 
speeds had lower kernel phase scatter. 


6. BINARY FITTING 

To search for companions in our data, we fit binary 
models to our kernel phases. Given the angular resolu¬ 
tion of our observations, a binary can be approximated 
as two delta functions, the Fourier transform of which 
is an analytic function. Equation |5]) gives the complex 
visibility for a binary with a companion separation s, po¬ 
sition angle P.A. (measured E of N), and brightness ratio 
b (in units of the primary’s brightness). 


V(u, v ) 



^ ^ gi-s(u-sin(PA)-\-v-cos(PA)) 


) 


( 5 ) 


The phase measured for a binary by a baseline with co¬ 
ordinates (u,v), is then the angle of Equation [ 5 ] 


4>{u, v) = tan 1 


/ Im(V(u,v))\ 
\Re(V(u,v))) ' 


( 6 ) 


To create a model set of kernel phases, we calculate the 
closure phases for each triangle using the locations and 
sizes of the mask sub-apertures, and sky rotation angles 
included in our observations. We then project the closure 
phases into kernel phases in the same way as we have 
done for the data. Due to the symmetry of the Fourier 
transform, each closure phase can correspond to one of 
two closing triangles. To keep the sign of our closure 
phases (and thus kernel phases) consistent between data 
and model, we sample the same closing triangles in both. 
Additionally, where necessary, we use observations of a 
known binary to calibrate the orientation of our detector 
on the sky. For the NaCo data, as in Huelamo et al. 
2011, we use observations of the binary 9 Ori C, taken 
in 2010 April. 

We perform y 2 fitting of the binary models to each in- 
dividual dataset using b oth a grid and nested sampling 
(Sivia & Skilling 2006). Our grid spans a range of pa¬ 
rameters in binary position angle (P.A.), separation (s), 
and contrast in magnitudes (A). A can be related to b, 
the brightness ratio, by the following: 


A - -2.5 log 10 (b) (7) 

Table [3] lists the ranges and spacings for each parameter. 
For each set of parameters, we calculated model kernel 
phases and a y 2 statistic. We used the best grid fit as an 
input for our nested sampling algorithm. 

Nested sampling involves first filling the parameter 
space with a large number (in our case, 100) of points, 
and calculating a likelihood (exp(— y 2 /2)) for each point. 
We then replace the lowest likelihood point with a ran¬ 
dom member of the remaining 99, and evolve it using a 


Markov-Chain to a higher likelihood region of the param¬ 
eter space. We repeat this process until the 100 points 
satisfy a convergence criterion, in our case, the scatter in 
the ensemble must be a small fraction (~ 0.1%) of the 
mean. We give one of the 100 nested sampling points 
the best grid fit as an initial value, and assign random 
values to the remaining 99. This is not necessary for 
fits in which there is one clear minimum, but it can help 
in preventing the nested sampling fit from converging to 
local minima. 

We report parameter errors calculated from a y 2 in¬ 
terval. After finding the minimum y 2 using nested sam¬ 
pling, we then scale all y 2 values for a grid of parameters 
so that the reduced y 2 of the best-fit model is equal to 1. 
With the scaled set of yf. values, we find all grid points 


1992) to place a lcr 


within Ay 2 of 3.53 (Press et al. 
error on the fit parameters. 

Bootstrapping often gives the most conservativ e esti¬ 
mates of par ameter errors. However, as noted in [Press] 
et al. (1992), the results of data fitting in Fourier space 
reiy heavily on all grid points being present, and while 
the ( u , v) coverage in the NaCo and MagAO data sets is 
good, it is by no means complete. For this reason, boot¬ 
strapped data sets do not fairly represent the noise in 
the data. We confirmed this by bootstrapping Gaussian 
noise sampled at the same (it, v) points as each dataset; a 
given noise realization’s best fit contrast ratio was much 
higher than the bootstrapped distribution would suggest. 

7. NOISE SIMULATIONS 

For each dataset, we fit simulated kernel phases to 
quantify our type I (false-positive) and type II (false¬ 
negative) errors. The contrast ratio and separation pa¬ 
rameters in a binary model act much like the amplitude 
and frequency of a sine wave; they take on non-zero val¬ 
ues when fit to noise. In order to determine the sepa¬ 
rations and contrast ratios that could be caused by the 
noise levels in these data sets, we simulated Gaussian 
kernel phases from distributions fit to the data. These 
are shown in Figure [3] for 1/ and 2013 Ks data, and in 
Appendix [A| for 2010 and 2011 Ks and H band data. We 
use ( u , v ) coverage and sky rotation identical to each ob¬ 
servation. We fit binary models to 1000 of these noise 
realizations and create a probability distribution from 
the best fits. Since the best fit position angle for the 
noise simulations is uniformly distributed, we create two- 
dimensional confidence intervals from the best fit separa¬ 
tions and contrasts. We then compare the best fit from 
our data to these confidence intervals. For example, if 
our best fit lies just outside the contour enclosing 95% of 
the best fits, there is a 5% chance that the fit is drawn 
from the distribution and thus a 5% chance that the fit 
resulted from noise alone. 

A second way of quantifying type I errors is to calcu¬ 
late the F statistic, the best fit y 2 divided by the null 
model y 2 , for each noise realization. Comparing the F 
statistic for a dataset’s best fit to a distribution of these 
simulated noise F statistics yields the probability that 
the best fit incorrectly rejects the null hypothesis. For 
example, if 95% of the simulated noise realizations have 
lower F statistics than a dataset’s best fit, there is a 95% 
probability that the data contai n no real signal. This 
procedure is outlined in detail in Protassov & van Dyk I 
( 2002 ). 
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TABLE 4 

Binary Fit Results 


Data Set 

P.A. (deg) 

Separation (mas) 

A L' 

P(FA) a (%) 

P(FA) b (%) 


& Hu,allowed 

2010 March 17: VLT/NaCo L' 






At" 1 

83 ±1 

88 ±§§ 

5-5 ±2:1 

11 

< 0.1 

2.63 

i 

At" 10 

82 ±f 

78 ±i 

5.3 ±0;| 

4 

< 0.1 

2.56 

i 

2011 March 14: VLT/NaCo L' 






At" 1 

92 ±H 

87 ±g? 

5.3 ±2 3 

25 

18 

2.54 

i 

At" 10 

92 ± 117 

94 ±lf 

5.4 ±°;f 

18 

36 

1.25 

i 

2012 March 8: VLT/NaCo L' 






At" 1 

-40 ±| 

32 ±| 8 

3.4 ± 1 ;i 

1.3 

< 0.1 

23.10 

> 4 

At" 10 

-38 ±| 

29 ±f 4 

0.4 ±g;l 

4.3 

< 0.1 

20.43 

4 

2013 March 25: VLT/NaCo L' 






At" 1 

83 ±f 

55 ±|1 

5-2 ±2:9 

0.9 

< 0.1 

25.82 

> 4 

At" 10 

82 ±i 

25 ±g° 01 

2.8 ±y 

3.4 

< 0.1 

17.35 

4 

2013 March 26: VLT/NaCo Ks 






At" 1 

74 ^0.3 

42 ±1 2 

5-2 ±?;i 

1.3 

< 0.1 

21.47 

4 

At” 10 

77 ±3 

51 ±® x 

s- 4 ±g;i 

0.4 

< 0.1 

10.64 

2 

2013 April 5: Magellan/MagAO/Clio2 L' 





At" 1 

112 ±1™ 

337 ± 1 tf 

5.8 ±i;| 

32 

2 

15.97 

4 

At- 10 

-131 ±1 

315 ±26 

5.8 ±1:| 

32 

4.5 

12.85 

3 


a Using distribution of noise simulation best fits 
b Using distribution of noise simulation F statistics 


To estimate our type II errors - the likelihood that 
the best fit to our data would be a false negative - we 
simulated observations with noise plus the signal from 
a companion. We generate 1000 Gaussian noise realiza¬ 
tions and fit a binary model to the simulated data. We 
take the fraction of fits where the input signal was not 
recovered to be our type II error. We also calculate F 
statistics from the noise + signal realizations, for com¬ 
parison with both the noise F statistics and the data’s 
best fit F statistics. 

Figure 4] shows the noise simulations for all 1/ and 
2013 Ks data sets with a At -1 calibration. The scattered 
points show the results of 1000 fits to Gaussian noise 
comparable to the scatter in each dataset, and the points 
with error bars show the best fit to each dataset. The 
color scale shows a probability distribution interpolated 
from the results of the simulations, and the contour lines 
indicate la, 2a, and 3cr confidence intervals. 

Figure [5] shows the false alarm probabilities for the 
same data sets in Figure [4] In each panel, the black line 
is a histogram of all of the F statistics for 1000 Gaussian 
noise realizations. The red, vertical lines indicate the 


F statistic for each dataset’s best fit; the intersection of 
these lines with the black histograms gives the probabil¬ 
ity that the best fit to the data resulted from noise. The 
green, cumulative histogram shows the distribution of F 
statistics for the 1000 noise + signal realizations. The 
intersection of the red line with this distribution gives 
the fraction of noise + signal realizations with fits that 
look less significant (higher F statistics) than our best fit. 
For Figures [4] and [5j the noise properties of the At -10 
calibration yielded comparable results (see Table [4]). 

While using a y 2 statistic to assign significance de¬ 
pends on the data error bars, which could be difficult 
to estimate properly, these simulations take into account 
only the scatter in the kernel phases. For this reason, in 
Section 8, we use the results of the noise simulations to 
estimate the significance of the best fit binary results. 

8. RESULTS 

Table [4] shows our results in chronological order. The 
first three columns list the results of binary fits to the ker¬ 
nel phases. Following the binary fits, the next columns 
list two false-alarm probabilities, the first calculated us- 
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NaCo L' 2010 
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NaCo L' 2013 



0 100 200 300 400 500 600 700 800 


NaCo L' 2011 



MagAO/Clio2 L' 2013 



Separation (mas) 


NaCo L' 2012 



NaCo Ks 2013 



Fig. 4.— Noise simulations for all L/ and NaCo 2013 Ks data sets. The scattered points show the best fits to 1000 noise realizations for 
each dataset, drawn from the Gaussian distributions shown in Figure [3] The color scale shows the probability distribution interpolated 
from the best fits, while the contours indicate lcr, 2cr, and 3cr confidence intervals. The bold point with error bars represents the best fit 
for each epoch. 


NaCo L' 2010 NaCo L' 2011 NaCo L' 2012 




NaCo L' 2013 



MagAO/Clio2 L' 2013 


NaCo Ks 2013 




F Xbest/~X.nidl 


Fig. 5.— False alarm testing for all L/ data sets and NaCo 2013 Ks observations. In each panel, the black, cumulative histogram shows 
the F statistics (best fit x 2 divided by null model % 2 ) for 1000 noise simulations. The intersection of the red, vertical line with the black 
histogram yields a false alarm probability. The green, cumulative histogram shows the F statistics from the 1000 noise + signal simulations 
carried out for each set of observations. The intersection of the red line with this histogram gives the fraction of noise + signal simulations 
that look less significant (higher F statistic) than the best fit. 
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Scan Index 


Fig. 6 .— Separation - contrast degeneracy. The subplots each 
show the evolution of one simulated kernel phase during the Ma- 
gAO/Clio2 observations. The different lines are binaries with iden¬ 
tical position angles, but different separations and contrast ratios. 
The blue solid li ne shows the signal fro m a companion similar to 
that published in Huelamo et al. |(2011) (s = 62 mas, A L' = 5.1), 
the red dashed line - shows s = 51 mas and A L' = 2.9, and the 
green dotted line shows s = 93 mas and A L' = 6.0. While slight 
differences exist between the binary models, they are not detectable 
given the accuracy of the data. 

ing the distribution of best fits to noise realizations and 
the second using the distribution of best fit F statistics 
(see Section [t|. The last two columns list the Ay 2 cor¬ 
responding to the Huelamo et al. (2011) binary for that 
dataset, alongside the corresponding confidence interval 
at which it is allowed. The bold values are binary fit re¬ 
sults for the At -1 calibrated data. We include the best 
fits from this calibration strategy in Section 9. 


8 . 1 . 2010 VLT/NaCo L' Data - Initial Detection 
For these data, we first fit our neighbor-calibrated 


strategy for the closure phases in 

Huelamo et al. 

( 2011 ) 

is most similar to this approach 

k see also 

Lacour et al. 


2011). We find a best fit binary with position angle of 
82 ± 7 0 , separation of 78 ± 2 § mas, and A V = 5.3 ± 
mag. This is comparable to the Huelamo et al. (2011) 
best fit - position angle of 78±1°, separation of 62±7 
mas (6.7 AU at 108 pc), and A I! of 5.1± 0.2 mag. Our 
fit param eter errors are substantia lly larger than those 
quoted in Huelamo et al. (20111, which were derived 


using \ 2 intervals from a binary fit to closure phases (see 


1 . : 
uld 


also Lacour et al. 2011). Closure phases have correlated 


errors, which could bias the parameter errors derived us¬ 
ing a simple \ 2 fit- Parameter errors also depend on 
whether the data error bars have been scaled such that 
the reduced \ 2 is equal to 1 , which is appropriate if you 
assume the model to be correct, and that underestimated 
error bars are causing a high reduced y 2 . Our large sep¬ 
aration and contrast error bars are consistent with the 
severe degeneracy between these two parameters. Figure 

f illustrates this degeneracy, with kernel phases plotted 
ir separations 0.5x, l x, and 1.5x the separation of the 
|Huelamo et al. ([2011 ) binary, varying the contrast ratio. 
The three models are nearly indistinguishable. For this 
reason, we believe our fit errors to be a more realistic 
representation of the parameter constraints. 

After the nearest-neighbor fit, we then fit the A t" 1 
calibrated kernel phases, finding a best fit with position 
angle of 83 ± g°, separation of 88± 2 g mas, and A L' = 
5 . 5 ± 2 '| mag. The A t," 1 calibration method resulted in 
kernel phases with lower scatter than the At -10 method 
- 0.63° versus 0.66°. Using a y 2 interval, the Huelamo et 
al. ( 2011 ) binary model is within ler of the best fit for 
both calibration methods. Appendix |B| shows our At" 1 
calibrated kernel phases with both the Huelamo et al. 
(2011) model and the best fit model from this work. The 
two models are nearly indistinguishable. 

Figure [4] shows that there is a small probability that 
the binary detection resulted from a random noise fluc¬ 
tuation. The best fit falls on a contour that encloses 89% 
of the simulation results; there is an 11 % chance that the 
2010 fit was the result of noise. The 2010 L' best fit is 
then significant at roughly the 2<r level. The distribution 
of F statistics, shown in Figure [5j indicates that the false 
alarm probability for this dataset is < 0 . 1 %, giving the 
fit ~ 3 <t significance. 


8 . 2 . VLT/NaCo 2011 II Data 

For these data, the At" 1 best fit has a position angle 
of 92±74°, a separation of 87±g7 mas, and a contrast of 
AL' = 5 . 3±2 3 mag. Compared to the At" 10 calibration 
(listed in Table [4|, the At" 1 weighting reduced the num¬ 
ber of outliers in the calibrated kernel phases, providing 
a tighter constraint on the binary fit parameters. Again, 
for both of these calibration methods, the Huelamo et al. 
(2011) model is within ler of the best fit. Appendix |B| 
shows our calibrated kernel phases with both the Hue¬ 
lamo et al. ( 2011 ) model and the best fit model from this 
work. 

Figure [4] shows the results of the noise simulations for 
this dataset. The distribution of best fits suggests that 
the fit to the data is significant at the ler level; the point 
with error bars falls on a contour which encloses 75% 
of the fits to noise, giving a 25% false alarm probabil¬ 
ity. This agrees roughly with the F statistic distribution, 
shown in Figure [5] The best fit F statistic, F = 0.900, 
gives an 18% probability of false alarm. This corresponds 
to a 1 — 2 cr confidence level. Here, the best fit is consistent 
with the distribution of noise + signal F statistics, which 
overlaps with the noise-only F statistics. The overlap of 
these two distributions indicates that, with the proper¬ 
ties of the NaCo 2011 data, noise alone can produce best 
fits that appear as significant as noise plus the Huelamo 
binary model. 
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8.3. VLT/NaCo 2012 V Data 

The A t~ 1 fit resulted in a position angle of —40±| o , 
separation of 32± 28 mas, and a contrast of 3.4±3'® mag. 
The kernel phases, with this work’s best fit and the Hue- 
lamo et al. (2011) companion signal over-plotted, are 
shown in Appendix [b] Using a y 2 interval, these data 
rule out the presence of the Huelamo signal at greater 
than 4cr. 

The results of noise simulations for the 2012 L' data 
are shown in Figure 0] The best fit lies on a contour 
that encloses 98.7% of the fits to noise, giving a 1.3% 
false alarm probability. This suggests that the best fit 
is significant at nearly 3er. The F statistics (see Figure 
[5]) also indicate that the fit is ~ 3 ct significant; less than 
TT.1% of F statistics are lower than that for the best fit 
(F = 0.579). 


8.4. VLT/NaCo 2013 L' Data 

The A t~ 2 best fit has a position angle of 83±f°, sep¬ 
aration of 55±25 mas, and a contrast of 5.2±5'g mag. 
Appendix [5] shows the kernel phases with our best fit 
model (realine) plotted alongside the Huelamo et al. 
(2011) binary model. The Huelamo et al. (2011) model 
is ruled out at ~ 4 sigma using a y 2 interval; the sepa¬ 
ration and contrast are within la of the Huelamo et al. 
(2011) binary, but the position angle is greater by 5.46°. 
However, adopting our larger parameter error bars, (see 
§8.1), the two models are consistent at la. 

Our best fit to these data lies on a contour enclosing 
99.1% of the points, giving a 0.9% chance that the best 
fit resulted from noise alone (see Figure |4|. This sug¬ 
gests that the best fit is significant at nearly 3cr. The 
F statistics, shown in Figure [5]) give a lower false alarm 
probability than this. The best fit F statistic (F = 0.586) 
is lower than all 1000 noise simulation F statistics, giving 
a < 0.1% false alarm probability. 

8.5. VLT/NaCo 2013 Ks Data 

The At -1 best fit parameters are position angle of 
74.45±g;g°, separation of 42.07±22 9 o7 mas j and contrast 
of 5.17±°;g7 mag. These data (scattered points) with our 
best fit model and the Huelamo et al. (2011) model are 
shown in Appendix[B] The error bars for this work’s best 
fit and the Huelamo et al. (2011) best fit overlap for both 
position angle and contrast. However, the best fit sepa¬ 
ration is smaller than the Huelamo et al. (2011) model. 
Given the large error bars, this difference is significant 
at less than 2cr. 

The best fit to these data falls on a contour that en¬ 
closes 98.7% of the noise simulation results (see Figure[4). 
This gives a 1.3% chance that the best fit was the result 
of noise. The best fit F statistic (F = 0.560) was lower 
than all of the simulated F statistics, shown in Figure [5j 
This gives a < 0.1% probability of the best fit resulting 
from noise. These estimates suggest that the fit to these 
data is significant at the ~ 3 a level. 

8 .6. Magellan/MagAO/Clio2 2013 Data 

The At -1 calibration for these data has a best fit with 
a position angle of 112±|g 6 °, separation of 337±(4 8 mas, 
and A L' = 5.8±o'!j mag. Both this and the At” 10 best fit 
(see Table]!]) are unreasonable in that binaries with these 


parameters should have been, but were not detected in 
the 2013 MagAO/Clio2 direct imaging data. 

Shown in Figure]!] the 2013 MagAO/Clio2 best fit lies 
on a contour that encloses 68% of the simulated fits, giv¬ 
ing a 32% probability that it could have resulted from 
noise - a ler result. For these data, this method does 
not agree well with the F statistic false-alarm estimation. 
The best fit F statistic, shown in Figure[5j is greater than 
only 2% of the simulated F statistics, suggesting that 
there is a 2% probability that the best fit was caused 
by noise. We speculate that this discrepancy could be 
caused by outliers in the data themselves. The null 
model x 2 for a set of kernel phases with non-Gaussian 
outliers will be greater than the null model for kernel 
phases drawn from a Gaussian distribution, which could 
reduce the F statistic. Figure [3] shows that the MagAO 
data could indeed have both outliers and a small non¬ 
zero mean kernel phase. These could both inflate the 
null model y 2 compared to that for Gaussian data, mak¬ 
ing the F statistic method less reliable as a false alarm 
probability estimator. 

Appendix [B] shows our calibrated kernel phases with 
the NaCo 2013 L 7 best fit and the best fit model. The 
NaCo 2013 L 7 best fit is allowed at 3cr by a secondary y 2 
minimum MagAO/Clio2 data. We estimated our type II 
errors following Section [TJ adding the 2013 NaCo V best 
fit to 1000 noise realizations. The results show that the 
probability of missing this signal, had it been present in 
these data, is 49.9%. Thus, the 2013 NaCo L 7 best fit is 
allowed by the MagAO/Clio2 observations. 

9. DISCUSSION 

The NaCo 2012 and MagAO/Clio2 2013 L 7 data sets 
have best fit binaries that are inconsistent with both the 
Huelamo et al. (2011) model and our best fit to the NaCo 
2013 L 7 data. However, simulations show that there is 
a non-zero chance that these fits resulted from noise - 
1.3% for the NaCo and 32% for the MagAO/Clio2 ob¬ 
servations. We simulated noise realizations (see Section 
[7]) to estimate our type II errors for both of these data 
sets, using the NaCo 2013 L 7 signal as the input binary 
model. The chance that we would have missed the binary 
signal in these data is 15.8% for NaCo 2012 and 49.9% 
for MagAO/Clio2 2013. 

Keeping these type II errors in mind, in the subsections 
concerning orbital motion and forward scattering we first 
assume that the tentative NaCo 2012 detection is reliable 
but that the MagAO/Clio2 non-detection is due to noise. 
We then discuss the results assuming that noise fluctua¬ 
tions led to a false detection in the NaCo 2012 data, while 
a signal compatible with the other NaCo data sets was 
actually present beneath the noise. We also consider the 
possibility that asymmetries caused by planet-disk inter¬ 
actions could have caused the observed kernel phases. In 
the last subsections, we discuss whether a chance align¬ 
ment or systematic error could masquerade as a compan¬ 
ion in the data sets where we detect a significant signal. 


9.1. Companion Orbital Motion 

Detecting orbital motion of a binary signal in multi¬ 
epoch data sets would confirm the claime d companio n 


from Huelam o et al. (2011). For example, Kraus & Ire¬ 
land (2012) detected a planet candidate in the LkCal5 
transition disk. Multi-epoch observations of this object 











New Spatially Resolved Observations of the T Cha Transition Disk 


13 



Fig. 7.— x 2 slices at the fixed contrast ratio for all L' observations as well as 2013 Na Co Ks observations, with filled contours at 1 to 
> 4cr confidence limits. The line indicates a circular orbit in the plane of the outer disk (see|01ofsson et al. |2013|). The xs show the initial 
position of the putative companion from our re-reduction of the 2010 NaCo L/ data, while the circles show the predicted position(s) of a 
planet on the orbit. We plot two since the planet could be orbiting in either direction. 


Example Orbits 



x NaCo 2010 

X NaCo 2011 

X NaCo + MagAO/Clio2 2013 



100 50 0 50 100 

a (mas) 

Fig. 8. — Example orbits consistent with 2010, 2011, and 2013 
NaCo best fits. The shaded regions show the lcr confidence in¬ 
tervals at the best-fit contrast ratio for the 2010 NaCo L/ (red), 
2011 NaCo L ' (blue), and 2013 NaCo L/+Ks (green) data sets. 
Colored points mark the predicted orbital positions at the times 
of the observations. Black ellipses mark the full orbits. The thick, 
gray line indicates the inner edge of the outer disk, as determined 
by Olofsson et al. (2013). 


revealed orbi tal motion of the compa nion at the level of 
~ 4° a year (Ireland & Kraus 2014), strengthening the 


case for the protoplanet explanation of the phase signal. 

We compare the predicted position of an orbiting com¬ 
panion to our fit results. Our best fit to the archival dis¬ 
covery data has a separ ation of 88 mas , which at T Cha’s 
distance of 108 ± 9 pc (Torres et al. ||2008| corresponds 
to 9.5 AU. We can use the orientation of the outer disk 
as well as the best-fit separation to predict the position 
of the planet in followup data sets. 

We first take a circular orbit at the same inclination 
(i = 58°) and po sition angle (—70°) as t he outer disk 


as determined by Olofsson et al. (2013). We choose 


the 2010 companion location as the initial orbital posi¬ 
tion, and then predict the projected separations for the 
followup data sets. Each panel in Figure [7] shows the or¬ 
bital prediction for a single observation over a y 2 slice at 
its best fit contrast ratio. While the predicted position 
is within lcr of the NaCo 2011 L' best fit, all other data 
sets, even those with lower false alarm probabilities such 
as 2013 NaCo 1/ and Ks, rule out the presence of a com¬ 
panion on this orbit. Thus, this scenario cannot explain 
the observations. 

If we assume the NaCo 2012 L ' signal is a real detec¬ 
tion, then we need to match its fitted position along with 
the fits from other epochs. This would require at least 
one full orbit to be completed between 2010 and 2013. 
However, the observed projected separations (s) at each 
epoch place a lower bound on the apocenter distance: 


a(l + e) > s. 


( 8 ) 
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Using this constraint and the fact that T Cha has a mass 
of 1.5Mq, we can check whether any orbits with a period 
of ~ 3 yr could produce the observed separations. Even 
as e approaches 1, an orbit around a 1.5M 0 star cannot 
have a period of less than 8.5 yr. This rules out orbital 
motion as the cause of the different position angle in the 

2012 L' data. 

If we assume that the 2012 L' NaCo data missed the 
signal found in the other data sets, we can ask what 
orbits would cause the best fits to the remaining obser¬ 
vations, which show motion compared to the Huelamo et 
al. (2011) model. For a grid of orbits, we calculated the 
predicted positions for the times of the 2010, 2011, and 

2013 observations. We compared these positions to the 
ler confidence intervals in position angle and separation 
from the binary fitting. We take the results of a simul¬ 
taneous fit to the 2013 NaCo and MagAO data as our 
constraint on the 2013 position. We find that inclined 
orbits, some of which cross into the outer disk, can pro¬ 
duce the observations. Orbits with the same inclination 
as the outer disk require eccentricities higher than 0.9 to 
reproduce the observations. Figure [8] shows two example 
orbits over the la confidence regions at the best fit con¬ 
trast for the three epochs. While a low-eccentricity orbit 
in the plane of the disk is inconsistent with the obser¬ 
vations, an orbit that is highly eccentric or substantially 
misaligned with the disk is compatible with the data. 
Since one would expect a young planet that is still ac¬ 
creting from the disk to be on a low-eccentricity, aligned 
orbit, these models do not seem physically likely. 


9.2. Forward Scattering from the Disk 


Olofsson et al. (2013) found that scattering by dust 
in the upper layers of T Cha’s outer disk could fit the 
observed phase signal nearly as well as a companion. In¬ 
deed, observations of another transition disk, FL Cha, 
showed that low mass companion and disk-scattering 
models could both explain th e closure phase sign al in 


this single epoch observation (Cieza et al. 2013). In 


contrast to the companion model, a constant scattering 
model would lead to a signal that does not vary in time. 

The NaCo 2011 and 2013 data support the hypothe¬ 
sis that a constant level of forward scattering could be 
responsible for the observations. The error bars for the 
best fits to these data sets are large enough that they 
overlap with those for our fit to the NaCo 2010 data. 
However, the best fit to the 2012 NaCo data is inconsis¬ 
tent with the NaCo 2010 best fit. If we assume that the 
2012 NaCo 1/ detection is not caused by noise, then vari¬ 
ability in the amount of scattered light would be required 
to explain the results. 

The intensity of scattered light in a protoplanetary disk 
is proportional to the luminosity of the star (e.g., Dong 
et al. | 2012} In oue et al. ||2008) |Dullemond fe INatta 
2003[) . Stellar variability would lead to variability in the 
intensity of the scattered light from the disk. However, 
this would not cause the A Li of a best fit companion to 
increase, since the ratio of the scattered light intensity 
to the stellar luminosity would remain constant. Addi¬ 
tionally, while T Cha is known to be quite variable in V 
band, due to changing extinction by circumstellar ma¬ 
terial, analysis of Spit zer spectra along with mid- and 
far-IR photometry by Schisano et al. (2009) indicates 
that it is not variable in the infrared. 


Changing the size distributions, cross sections, or mass 
fractions of the various dust grain species could change 
the amount of f orward scattering relative to the stellar 
bright ness (e.g., D’ Alessio et al. 1998 Pollack et al. 


1985). Changes to the geometry of the disk, such as 
the height of the outer disk’s in ner wall, coul d also al¬ 


ter th e scattering intensity (e.g., Dullemond & Dominik 
2004). Some young stars exhibit variability in scattered 
light due to geometric changes in the in ner disk that then 
shadow the outer parts of the disk ([Wisniewski et al. 


2008] Sitko et al. 2008 Bans & Kmgl |20I2 e.g.,). These 
changes can occur on the timescale of weeks to months 
(the dynamical timescale of the inner disk). However, 
geometric changes at the radius of the outer disk would 
take place on the timescale of multiple years. The vis¬ 
cous timescale, an estimate of the time it takes for disk 
material to shear out, at 9.5 AU, is greater than 100 yr 
for reasonable valu es of the viscosity parameter (Shakura 
& Sunyaev ] 1973) and scale-height to radius ratio? The 


dynamical timescale (D^ 1 ) at the inner edge of the outer 
disk (12 AU) is 5.41 yr. Thus, disk geometry changes 
at the radius of the outer disk cannot explain the 2012 
NaCo L' best fit. 

Scattered light from the upper layers of the outer disk 
should be brighter at Ks than at L'. This is due to the 
fact that larger grains will settle toward the disk mid¬ 
plane and thus contribute less to the to tal am ount of 
scatt ered light than smaller grains (e.g., |Natta et al. 
2007). Assuming that the detections in the JNaCo 2012 
and MagAO data are due to noise, and that forward 
scattering is causing the other observed signals, one may 
naively expect the 2013 kernel phase signal to be greater 
at Ks than at I/. In contrast, we observe similar signal 
amplitudes at both wavelengths. However the scattered 
light may arise from an extended region, significantly re¬ 
solved by our observations, in which case such a simple 
interpretation may not apply. To test this, we simu¬ 
lated scattered light images of T Cha’s outer disk usin g 
the radiative transfer code Hyperion (Robitaille |2011). 
We produced images comparable to those published in 
Olofsson et al. (2013), using the same dust properties 
and disk parameters as their best fit. For this model, 
the mean kernel phase signal was 0.29° at Ks and 0.27° 
at L'. Thus, having similar kernel phase signals for Ks 
and L' cannot rule out the scattering scenario. We also 
calculated % 2 values using the 2013 NaCo observations, 
and find that the reduced \ 2 for the scattering model is 
~ 5.1. The binary model gives a better fit to the data, 
but forward scattering can produce kernel phases similar 
to the observations. 


9.3. Optically Thin Disk Asymmetries 

While the scattering scenario leads to asymmetries in 
the outer disk, resulting in non-zero kernel phases with¬ 
out the pr esence of a com panion, hydrodynamic simula¬ 
tions (e.g., Fouchet et al. 2010) suggest that disk-planet 
interactions can cause asymmetric structures in the opti¬ 
cally thin dust within transition dis k cavities. Observa- 
tions in both the inf rared (e.g., Muto et al. |2012 ) and the 
sub millimeter (e.g., Isella et al. T30T5f]Ferez = efal. |2014|) 
have confirmed the presence of such asymmetries. Fur- 
thermore, recent NRM observations of the transition disk 
V1247 Orionis revealed phase signals whose best fit bi- 
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nary parameters changed significantly with wavelength. 
This indicated that the underlying structure was not a 
simple companion, but as ymmetric optica lly thin mate¬ 


rial within the disk gap (Kraus et al. 2013). Image 


reconstruction would allow us to look for asymmetries 
in T Cha’s disk gap. We will present image reconstruc¬ 
tion in a systematic way in a f uture paper, bu t a pre ¬ 


liminary effort (using MACIM; Ireland et al. (2006)) 


suggests a trefoil structure within the cleared region of 
the transition disk, perhaps compatible with disk-planet 
interaction models. 


9.4. Chance Alignment 

We investigate the probability that the chance align¬ 
ment of a foreground or background object would cause 
a companion signal in the 2010 NaCo dataset, and that 
its proper motion would be mistaken for orbital motion 
between 2010 and 2013. Since T Cha has a proper mo¬ 
tion of -39.61 mas yr -1 in right ascension and -9.87 mas 
yr -1 in declination, stationary objects in the foreground 
or background could appear to move like orbiting com¬ 
panions. 

The number of chance alignments in a field of view 
with area Apov is 


nalign — ApQyYj^ (9) 

where £ is the surface density of stars along the line of 
sight. We use the extent of T Cha’s disk gap to define our 
field of view. From Olofsson et al. 2013, the gap extends 
from 0.17 AU to 12 AU, subtending ~ 0.11” at 108 pc. 
To estimate the surface density of stars along the line of 
sight, we queried Two Micron All Sky Survey for all stars 
within 1 degree of T Cha, with AKs lower than 7 (Ks of ~ 
7—14 mag). This gives a stellar surface density of 11544 
stars deg -2 . The probability that a chance alignment 
with a foreground or background object would cause a 
companion signal in the NaCo 2010 data is then ~ 7.0 x 
10 -6 . 


9.5. Systematic Errors 

The 2010, 2011, and 2013 best fits have a low, 0.003%, 
probability of all resulting from noise, assuming the ob¬ 
servations are independent. In this section, we discuss 
the possibility that the signals could be caused by some 
systematic effect. 

Systematic errors could be possible since T Cha is such 
a southern target (5 ~ —79°), and thus transits at 35.27° 
at the VLT (cp = —24.63°) and slightly higher, 39.90° at 
Magellan (cp = —29.26°). While our calibrator tests (see 
Section [5]) showed that it is unlikely that the calibrators 
injected a signal into the T Cha kernel phases, this does 
not rule out all AO-related systematic effects. If we as¬ 
sume a reliable non-detection in the 2013 MagAO/Clio2 
data, this could suggest that a systematic error caused 
the same signal to be present in all of the NaCo data sets. 
This could be due to T Cha’s lower transit elevation as 
observed from the VLT. 


The similar position angles, but different separations 
of the L' and Ks band best fits (see Figure [7]) also sug¬ 
gest that systematic errors could be an issue. The best 
fit separation for the 2013 NaCo Ks observations is ~ 42 
mas, while the L' best fit separation is ~ 55 mas. The 
size of an observ ed speckle patter n is proportional to the 
wavelength (e.g., Marois et al.]2000), and thus we would 
expect the ratio of a speckle’s position in Ks to its posi¬ 
tion in L' to be ~ 0.57. The ratio of the NaCo 2013 Ks 
and L' separations is 0.77. However, the large separation 
error bars prevent us from placing tight constraints on 
the relative separations. 


10. CONCLUSIONS 

We presented multi-epoch observations of the T Cha 
transition disk, taken using VLT/NaCo and Magel- 
lan/MagAO/Clio2, in L', Ks, and H bands. Out of the 
nine data sets, three are too noisy to detect signals at the 
level of the companion candidate published in Huelamo 
et al. (2011); these are the 2010 and 2011 NaCo Ks and 
2013 NaCo H band data. 

We find companion parameters comparable to those 
published in Huelamo et al. (2011) in our re-reduction 
of the initial discovery data. Furthermore, we find com¬ 
parable binary parameters to the companion candidate’s 
in the 2011 NaCo L' dataset, and tentative evidence for 
radial motion by the time of the 2013 NaCo and Ma- 
gAO/Clio2 observations. Assuming that noise in the 
2012 NaCo and 2013 MagAO/Clio2 data led to non¬ 
detections of the signal (simulations show such false non¬ 
detections to be fairly probably in these data sets), highly 
eccentric or misaligned orbits could result in a signal con¬ 
sistent with the observations. 

Scattered light from T Cha’s outer disk could provide 
another explanation for the observed kernel phases, al¬ 
though a binary model gives a slightly better fit to the 
data. Preliminary image reconstructions also suggest 
an asymmetric structure, perhaps consistent with disk- 
planet interactions, as the source of the observed sig¬ 
nals. Lastly, the ratio of the NaCo 2013 Ks and L' best 
fit separations argues for the possibility that challenges 
associated with AO correcting a dim, southern source, 
could cause a systematic error that would masquerade 
as a close in companion. The detection of a secondary 
minimum in the MagAO/Clio2 data at the same position 
as the NaCo detections is encouraging, but follow-up ob¬ 
servations with higher signal to noise are required to rule 
out the possibility of systematic errors. 
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NaCo Ks 2010 



NaCo Ks 2011 



NaCo H 2013 



Fig. 9.— Sky rotation comparison in ( u,v ) space for 2010 and 2011 NaCo Ks data as well as 2013 NaCo H data. The NaCo 2010 and 
2011 Ks data had changes in sky rotation of ~ 19° and ~ 38°, respectively. The sky rotation in the NaCo 2013 H data was ~ 39°. 


Snapshot Kernel Phase Scatter 



Fig. 10.— Normalized histogram of uncalibrated kernel phases for 2010 and 2011 Ks as well as 2013 H band data. For a subset of each 
dither (taken so that equal amounts of integration came from each observation), we subtract the mean kernel phase from each individual 
measurement to generate the histogram shown above. The snapshot errors for these data sets are significantly larger than those for the 
other six. 


APPENDIX 

A. ADDITIONAL KS AND H BAND DATA SETS 


Figure [9] shows the sky rotation coverage for the 2010 and 2011 NaCo Ks and 2013 NaCo H band data sets. The 
comparison of the snapshot errors (see ® is displayed in Figure 10 The calibrated kernel phase histograms, with the 
Gaussian distributions used to generate the noise simulations discussed in Section [TJ are shown in Figure El 

A.l. 2010 VLT/NaCo Ks Data - Published Non-detection 

Table [5] lists the results of binary fits to these data as well as the 2011 NaCo Ks and 2013 NaCo H band data 
sets. The AK 1 calibration strategy produced a best fit position angle of 122±f3 7 °, separation of 268±|g 2 mas, and a 
contrast of AA's = 4.8 ±q'5 mag. Appendix |b| shows these data alongside this work’s best fit and the Huelamo et al. 
(2011) binary model. 

For these data, there is a very high probability that the best fit resulted from noise, and that a companion signal 
would not have been recovered. The best fit to these data lies on a contour that encloses 29% of the noise simulation 
results (see Figure 12). There is thus a 71% chance that the best fit resulted from noise. The F statistics agree with 
this result; 71.5% of the F statistics are lower than that measured for the best fit, F = 0.843 (see Figure [l3]) . Both of 
these suggest that the fit to this dataset is consistent with noise. 
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A t~ l Calibrated Kernel Phases 



Fig. 11.— Normalized histograms of calibrated NaCo 2010 and 2011 Ks as well as NaCo 2013 H kernel phases, with their best fit Gaussian 
distributions over plotted. These calibrated kernel phases have much higher scatter than the other data sets. The Gaussian distributions 
were used to generate noise realizations for the simulations described in Section [7] 


TABLE 5 

Binary Fit Results 


Dataset PA (deg) Separation (mas) A L' P(FA) 1 (%) P(FA) 2 (%) P(miss) (%) A x 2 Hu &Hu,allowed 


2010 July 1: VLT/NaCo Ks 


At' 1 

122 ±l% 7 

268 ± 412 

4 Q _|_1.3 

4,0 =^0.5 

71.0 

71.5 

93.6 

10.40 

3 

At- 10 

124 ±\H 

259 

4 o J-1.8 
^0.6 

36.0 

92.6 

95.5 

8.19 

3 

2011 March 15: VLT/NaCo Ks 







At' 1 

ryry _|_0.03 
± 165 

358 ±g 42 

4 o -l0.6 
=t 0.4 

68.0 

9.5 

83.5 

9.44 

3 

Ai- 10 

22 i 1 

359 ±1 

4 -i ±8:1 

31.7 

5.7 

82.9 

9.31 

3 

2013 March 27: VLT/NaCo H 







At' 1 

-27 ±? 

31 

3.8 ±°;i 

1 

< 0.1 

84.1 

32.01 

> 4 

Ai- 10 

-28 ±| 

29 ± 18 

3.6 ±8;§ 

0.8 

< 0.1 

98.1 

30.10 

> 4 


a Using distribution of noise simulation best fits 
k Using distribution of noise simulation F statistics 


For the 2010 and 2011 NaCo Ks and 2013 NaCo H band data sets, we estimate our type II errors following the 
procedure outlined in (j7j As the input signal, we add companions with the same separation (62 mas) and contrast 
(A = 5.1) as the Huelamo et al. (2011) binary. To account for biases due to sky rotation coverage, we vary the input 
position angle of the companions. Of the 1000 noise + signal realizations, 936 resulted in erroneous best fits. This 
gives a 93.6% chance that, with the noise present in the 2010 Ks data, we would not have recovered the companion 
signal had it been present. We thus cannot assign a non-detection to this dataset with a high degree of confidence. 

A.2. VLT/NaCo 2011 Ks Data 

The At' 1 calibration yielded a best fit position angle of 22±°g° 5 3 °, separation of 358±| 42 mas, and a contrast of 
AKs = 4.2 ±q;4 5 mag. See Tablets for At -10 best fit parameters. Appendix [ b] shows these data with this work’s best 
fit and the Huelamo et al. (2011)binary model plotted. 

Figure [l2| shows the results of noise simulations for the 2011 NaCo Ks data. The best fit falls on a contour that 
encloses 32% of the noise simulations. This gives a 68% chance that the best fit was caused by noise alone. The F 
statistics (see Figure [l3j), however, give a smaller false alarm probability; the best fit F statistic (F = 0.837) is higher 
than only 9.5% of the simulated F statistics. This gives a 9.5% chance that the best fit resulted from noise. 
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NaCo Ks 2010 NaCo Ks 2011 



Separation (mas) 


NaCo H 2013 

I Best Fit 



Fig. 12.— Noise simulations for all H and Ks data sets. Th e scattered points show the best fits to 1000 noise realizations for each dataset, 
drawn from the Gaussian distributions shown in Figure EH The color shows the probability distribution interpolated from the best fits, 
while the contours indicate la, 2 a, and 3 a confidence intervals. The bold point with error bars represents the best fit for each epoch. 


NaCo Ks 2010 



NaCo Ks 2011 



NaCo H 2013 



Fig. 13.— False alarm testing for 2010 and 2011 Ks and 2013 H band observations. In each panel, the black, cumulative histogram shows 
the F statistics (best fit y 2 divided by null model y 2 ) for 1000 noise simulations. The intersection of the red, vertical line with the black 
histogram yields a false alarm probability. The green, cumulative histogram shows the F statistics from the 1000 noise + signal simulations 
carried out for each set of observations. The intersection of the red line with this histogram gives the fraction of noise + signal simulations 
that look less significant (higher F statistic) than the best fit. 


The type II error estimation gives an 83.5% chance that the companion candidate signal would be lost under the 
noise. Of the 1000 realizations, 835 resulted in best fits that were inconsistent with the input signal. Thus, we cannot 
rule out the Huelamo et al. (2011) signal in these data. 


A.3. VLT/NaCo 2013 H Data 

The At -1 best fit to these data has a position angle of -27±y 0 , separation of 29±g 8 , and contrast of 3.6±o g mag. 
Table [ 5 ] lists the fit parameters for the At -10 calibration. See Appendix |b| for a plot of these data, along with this 
work’s Dest fit model and the Huelamo et al. (2011) companion model. 

Figure [12] shows the results of noise and noise + signal simulations for the 2013 H band data. The best fit to these 
data, compared to the distribution of fits to noise, suggests that there is only a 1% chance that the best fit is the 
result of noise. The F statistics also indicate that the false alarm probability is < 0.1% (see Figure 13); the best fit 
F statistic (F = 0.725) is lower than all simulated F statistics. However, the kernel phase histogram (see Figure 11) 
shows a significant number of outliers. The distribution is also not centered on zero. Both of these characteristics 
would cause the distribution of best fits and F statistics to underestimate the false alarm probability. 

Of the 1000 noise + signal simulations, we did not recover the input signal in 841 realizations. This gives an 84.1% 
chance that we would have missed the companion signal, had it been present under the noise. Thus, we cannot rule 
out the Huelamo et al. (2011) signal with confidence using the H band data. 
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Figures 14 - 22 show plots of kernel phase versus scan index, a proxy for sky rotation angle, for all observational 
epochs. In each plot, the grey points show the data and the red solid line marks this work’s best fit. Figures 14 - 19 
show data sets with low enough noise levels to include in our Discussion. Figure 19, which displays the MagAO/Clio2 
data, shows the NaCo 2013 L' best fit in green, since the two data sets were taken close in time to one another. All 
other figures in this section show the Huelamo et al. (2011) model in green. Figures 20 - 22 show data sets with high 
enough scatter to wash out signals of interest, and thus were not included in our Discussion. 


NaCo 2010 L' At 1 Calibration 



Scan Index 


Fig. 14.— Kernel phase data for 2010 NaCo L' observations, shown with best fits from this work (red line) and Huelamo (dashed green 
line). Each subplot corresponds to a single linear combination of closure phases plotted against scan index (a proxy for sky rotation angle). 
The error bars plotted are unsealed, while our parameter constraints are derived using error bars scaled such that the reduced y 2 of the 
best fit model is equal to 1. The Huelamo et al. (2011) model is allowed within ltr of this work’s best fit. 


C. KERNEL PHASE PROJECTION 


Our kernel phase projection is similar to that in 


the ways in which N a apertures (</>) are combined to yield iN f 

$ = A 


Martinaclie (2010). We begin with the matrix A, which describes 
phases (<h): 

<t>- (Cl) 


This equation can be modified for observations of a source with intrinsic signal, assuming that the source phase simply 
adds to the instrumental phase: 

$ = A-^ + $ 0 - (C2) 

In order to eliminate the instrumental phase, we are searching for a matrix, K, such that: 

K A = 0. (C3) 

We can use singular value decomposition to find K. We decompose A 1 in the following way: 

A T = U • W • V T (C4) 


where U is an N a x N p column-orthogonal matrix, W is an N p x N p diagonal matrix with either positive or zero 
elements, and V is an N p x N p orthogonal matrix. The columns of V corresponding to zero W-values are filled into 
the rows of K. We then build K from a matrix, T, which describes how to combine phases into closure phases (“hep). 

$ cp = T- A-0+T-$ o = T-$ 0 (C5) 


We find B, such that: 


B T = K 


(C6) 


and 


$ fc =B-T-A-0 + B- T-$o=K-$o 


(C7) 
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Fig. 15.— Kernel phase data for 2011 NaCo L/ observations, shown with best fits from this work (red line) and Huelamo 2010 (dashed 
green line). Each subplot corresponds to a single linear combination of closure phases plotted against scan index (a proxy for sky rotation 
angle). The error bars plotted are unsealed, while our parameter constraints are derived using error bars scaled such that the reduced x 2 
of the best fit model is equal to 1. The Huelamo et al. (2011) model is allowed within la of this work’s best fit. 
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Fig. 16. — Kernel phase data for 2012 NaCo L' observations, shown with best fits from this work (red line) and Huelamo (dashed green 
line). Each subplot corresponds to a single linear combination of closure phases plotted against scan index (a proxy for sky rotation angle). 
The error bars plotted are unsealed, while our parameter constraints are derived using error bars scaled such that the reduced y 2 of the 
best fit model is equal to 1. 

In the following equations, for any matrix M, M ~ ight and M ; ~j t represent the right and left generalized inverses, 
respectively. These are used to invert non-square matrices; a right inverse is required when a matrix has full row rank 
but does not have full column rank, and a left inverse when a matrix has full column rank but does not have full row 
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Fig. 17.— Kernel phase data for 2013 NaCo L/ observations, shown with best fits from this work (red line) and Huelamo (dashed green 
line). Each subplot corresponds to a single linear combination of closure phases plotted against scan index (a proxy for sky rotation angle). 
The error bars plotted are unsealed, while our parameter constraints are derived using error bars scaled such that the reduced % 2 of the 
best fit model is equal to 1. 
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Fig. 18. — Kernel phase data for 2013 NaCo Ks observations, shown with best fits from this work (red line) and Huelamo (dashed green 
line). Each subplot corresponds to a single linear combination of closure phases plotted against scan index (a proxy for sky rotation angle). 
The error bars plotted are unsealed, while our parameter constraints are derived using error bars scaled such that the reduced % 2 of the 
best fit model is equal to 1. 


rank. 


M ~? ght = M T • (M • M T ) -1 

(C8) 

= ( mT • M ) _1 • mT 

(C9) 


K has the dimensions N& x N p where is the number of (linearly independent) kernel phases and N p the number 
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Mag AO L' 2013 At 1 Calibration 



Fig. 19.— Kernel phase data for 2013 MagAO L/ observations, shown with best fits from this work (red line) and the NaCo 2013 L/ best 
fit (dashed green line). Each subplot corresponds to a single linear combination of closure phases plotted against scan index (a proxy for 
sky rotation angle). The error bars plotted are unsealed, while our parameter constraints are derived using error bars scaled such that the 
reduced % 2 of the best fit model is equal to 1. 
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Fig. 20.— Kernel phase data for 2010 NaCo Ks observations, shown with best fits from this work (red line) and Huelamo (dashed green 
line). Each subplot corresponds to a single linear combination of closure phases plotted against scan index (a proxy for sky rotation angle). 
The error bars plotted are unsealed, while our parameter constraints are derived using error bars scaled such that the reduced % 2 of the 
best fit model is equal to 1. 


of Fourier phases (N& < N p ). K has a right inverse. 


B T K 


right 


= K ■ K 


right 


= i. 


(CIO) 


















































New Spatially Resolved Observations of the T Cha Transition Disk 

NaCo 2011 Ks At -1 Calibration 


23 



Fig. 21. — Kernel phase data for 2011 NaCo Ks observations, shown with best fits from this work (red line) and Huelamo (dashed green 
line). Each subplot corresponds to a single linear combination of closure phases plotted against scan index (a proxy for sky rotation angle). 
The error bars plotted are unsealed, while our parameter constraints are derived using error bars scaled such that the reduced % 2 of the 
best fit model is equal to 1. 


NaCo 2013 H At 1 Calibration 


Oh 



Scan Index 


Fig. 22. — Kernel phase data for 2013 NaCo H observations, shown with best fits from this work (red line) and Huelamo (dashed green 
line). Each subplot corresponds to a single linear combination of closure phases plotted against scan index (a proxy for sky rotation angle). 
The error bars plotted are unsealed, while our parameter constraints are derived using error bars scaled such that the reduced % 2 of the 
best fit model is equal to 1. 


Since B has the dimensions N& x N cp , 


(Nfc < N cp ), B has a right inverse, which, according to CIO 


is the following: 


B 


-l 

right 


= T K 


-l 

right 


(Cll) 


B B 


-1 

right 


= i 


and 


(C12) 







































































24 


Sallum et al. 


Furthermore, the left inverse of B ri J u can be calculated to find the B in Equation (C12). 
So: 

g_ i \—l _ /it' t/— i \—l 


(B 


rightheft ’ ^right) left 


(CIS) 


This satisfies both Equations (C6| and (C7|. We find the kernel phase covariance matrix, Cfc, using the closure phase 
covariance matrix C cp in the following way: 

C fc = B • C cp ■ B t . (C14) 

We can take our kernel phase variances to be the diagonal entries of Cfc. 
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